Structure, reactivity and electronic properties of V-doped Co clusters 
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Structures, physical and chemical properties of V doped C013 clusters have been studied in detail 
using density functional theory based first-principles method. We have found anomalous variation in 
stability of the doped clusters with increasing V concentration, which has been nicely demonstrated 
in terms of energetics and electronic properties of the clusters. Our study explains the nonmonotonic 
variation in reactivity of Coi3- m V m clusters towards H2 molecules as reported experimentally [J. 
Phys. Chem. 94, 2744 (1990)]. Moreover, it provides useful insight into the cluster geometry and 
chemically active sites on the cluster surface, which can help to design better catalytic processes. 



I. INTRODUCTION 

The interest in the studies of clusters is largely because 
of their possible technological applications which include 
the possibilities of developing novel cluster-based materi- 
als using the size dependence of their properties. Doping 
of clusters is an important possibility in this direction. In 
recent times the fabrication of alloy clusters of different 
sizes with well denned, controlled properties by varying 
the composition and atomic ordering, has caught consid- 
erable attention!^ Bimetallic alloy clusters have been 
known and exploited for last few years in various ap- 
plications, specially in the catalytic reactions.^ Varying 
the ratio of the two constituents, one can alter the surface 
structures, compositions and segregation properties!^ In 
this way, it is possible to tune the chemical reactivity at 
the surface of an alloyed cluster Considering the huge 
possibility of using nano clusters in catalysis, the study 
of cluster reactivity has become a n interdisciplinary topic 
of present day research! 9 * 10 * 11 * 12 ^ 

About two decades ago, Nonose et al measured the re- 
activity towards H2 of bimetallic Co n V m (n > m) clusters 
using laser vaporization technique and reported strong 
cluster size and composition dependence.^ Both V and 
Co are 3d transition metals. The substitution of Co by V 
atoms, one by one, should increase the reactivity of the 
alloy cluster towards H2 molecules, as V, an early transi- 
tion metal, has a high reactivity towards H2 in contrast 
to Co which has relatively low reactivity. 14 Interestingly, 
it was found that the reactivity increased gradually, as 
expected, with the substitution of Co atom by V atom 
one by one in the Co n clusters (n < 13), but for the C013, 
there was a remarkable decrease in reactivity when a sin- 
gle Co atom was substituted by a V atom. However, the 
reactivity increased as the number of exchange V-atoms 
increased further up to m = 3, while the fourth V-atom 
substitution did not increase the reactivity any more. In 
view of the high reactivity of elemental V, this sudden 
drop in reactivity of the C012V cluster was rather sur- 
prising. The authors speculated a plausible icosahedral 
structure for the C012V cluster with the active V atom 



at the cage center. Therefore, the V atom, being shielded 
geometrically from H2 by the twelve surface Co atoms, 
might have less chance to interact with H-atoms. The 
chemisorption reactivity of cationic Coi3_ m V+ clusters^ 
and anionic Coi3_ m V~ cluster^ also shows similar type 
of variation as that of neutral Coi3_ m V m clusters, which 
hints towards the dominant effect of geometric structure 
as compared to the electronic structure. For clusters, 
the ionization potential (IP) depends on the position of 
the highest occupied molecular orbital (HOMO) level. 
For the pure metal clusters, Whetten et al have postu- 
lated that the reaction rate for cluster-H2 dissociative 
chemisorption is determined by the charge transfer from 
the HOMO to the lowest unoccupied molecular orbital 
(LUMO) of the reactant gas H2, and an anti-correlation 
between IP and reaction coefficient could be observed.^ 
That means lower value of IP corresponds to higher re- 
activity and vice versa. However, the ionization ener- 
gies of the Co n V m clusters, measured by Hoshino et al 
using photo-ionization spectroscopy show no such anti- 
correlationj^ again demonstrating the importance of ge- 
ometrical structure. A rigorous first-principles study in 
terms of geometric and electronic effects is therefore very 
much needed to understand this anomalous nature of re- 
activity of Coi3_ m V m (m = 0-4) clusters. In this report, 
we have carried out an ab-initio theoretical study on V 
doped C013 clusters. The whole study can be divided 
into three major parts : first part consists of an exhaus- 
tive search for the minimum energy structures (MES) for 
cluster of each composition, followed by stability analy- 
sis of these MES in terms of various physical quantities 
in the second part, while the final part includes inves- 
tigation of chemisorbed structures and understanding of 
cluster reactivity. 

II. COMPUTATIONAL DETAILS 

The calculations have been performed using den- 
sity functional theory (DFT), within the pseu- 
dopotential plane wave method as implemented in 
VASP codeP^ We have used the projector aug- 
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mented wave (PAW) pseudopotentialj 20 * 21 * and the 
Perdew-Bruke-Ernzerhof (PBE) exchange-correlation 
functional 22 for spin-polarized generalized gradient ap- 
proximation (GGA). The 3d and 4s electrons were 
treated as the valence electrons for the transition metal 
elements and the wave functions were expanded in the 
plane wave basis set with the kinetic energy cut-off of 
335 eV. The convergence of the cluster energies with re- 
spect to the cut-off value has been checked. Recipro- 
cal space integrations were carried out at the T point. 
Symmetry unrestricted geometry and spin optimizations 
were performed using the conjugate gradient and the 
quasi-Newtonian methods until all the force components 
were less than a threshold value of 0.005 eV/A. In order 
to ensure the structural trends found in the optimized 
structures, we have also carried out Bohn-Openheimer 
molecular dynamics 2 ^ within LDA for few specific clus- 
ters. Simple cubic super-cells were used for cluster cal- 
culations with the periodic boundary conditions, where 
two neighboring clusters were kept separated by at least 
12 A vacuum space to make the interaction between the 
cluster images negligible. The cohesive energy (E c ) of a 
Co n V m alloy cluster is calculated with respect to the free 
atoms as 

£ c (Co n V m ) = mE(V) + n£(Co) - E(Co n V m ) (1) 

where E(Co n V m ), E(Co) and E(V) are the total energies 
of Co n V m cluster, an isolated Co atom and an isolated 
V atom respectively. One can also define the cohesive 
energy with respect to the Co and V bulk phases at equi- 
librium instead of isolated atoms. However, the cohesive 
energy variation with increasing V-concentration follows 
the same trend irrespective of the way of cohesive energy 
definition. The second difference in energy for fixed size 
(n + m = constant) and variable composition is defined 
as 

A 2 £(Co n V m ) = ^(Co n+1 V m _ 1 )+^(Co n _ 1 V m+ i) 
- 2£(Co n V m ) (2) 

It gives the relative stability of alloy clusters having 
nearby compositions. 

III. STRUCTURE 

Since chemical reaction with clusters takes place on its 
surface, the atomic arrangement and the composition of 
cluster surface play an important role in chemical reac- 
tivity. Therefore, the first step towards the theoretical 
modeling of clusters is to determine their ground state 
structure. In our earlier workpl we studied the struc- 
ture and magnetism of the pure Co clusters in detail. To 
obtain the MES having the optimized geometry as well 
as the optimized magnetic moment, we considered sev- 
eral probable starting geometries having closed packed 
atomic arrangement and allowed each of the geometries 
to relax for all possible collinear spin configurations of the 



atoms. We follow the same way of structural optimiza- 
tion here for the V doped Co clusters ( Coi3_ m V m ; m = 
1-4). However, the situation for alloy clusters is quite 
cumbersome as one has to deal with a large number of 
starting geometries because of the pr esence of homotops 
(the term was first coined by JellinelP^]) . The homo- 
tops have the same number of atoms, composition and 
geometrical structure, but differ in the arrangement of 
the doped atoms. As for example, a single geometrical 
structure of an A n B m alloy cluster with fixed number of 
atoms (TV = m + n) and composition (m/n ratio), will 
give, in principle, -^p-\ homotops. Some of them, how- 
ever, may be symmetry equivalent and the number of 
the inequivalent homotops will be somewhat less than the 
above mentioned value. So the variety of structures in al- 
loy clusters is much richer than that of the pure clusters 
and the potential energy surface of even a small cluster 
of few tens of atoms is of enormous complexity. 

It is generally found that the transition metal clusters 
prefer compact geometries to maximize the interaction 
between the rather localized d orbit alsPS There are three 
most compact and highly coordinated structures for 13 
atoms cluster: icosahedron, cub-octahedron and hexago- 
nal closed packed (HCP) geometries. We have therefore 
considered these three geometries as the most probable 
starting structures. For the pure C013 cluster, we found 
that the MES is a distorted HCP geometry with total 
magnetic moment of 25 hb and total cohesive energy of 
42.63 eV as mentioned in our previous work. 24 The struc- 
ture has 22 triangular faces and 33 edges (cf. Fig. [I]). 
This hexagonal growth of the pure Co clusters is quite 
unusual compared to the trend seen in the clusters of 
its 3d neighbors like Mn, Fe, N i, which preferably adapt 
an icosahedral pattern PHESM] Further study is needed to 
explore the u nusual structural pattern of small pure Co 
clusters!- 31 * 32 * Another distorted structure of HCP motif 
with total magnetic moment of 27 \±b and 0.14 eV above 
the minimum energy state, is the first isomer. The op- 
timal icosahedral structure of total spin 31 \ib (struc- 
ture having 20 triangular faces and 30 edges) is 0.17 eV 
higher than the minimum energy state and emerges as 
the second isomer. The third isomer is a distorted cub- 
octahedron with total magnetic moment 25 hb and it 
is 0.22 eV above the minimum energy state. The struc- 
tures of the ground state, second and third isomers are 
shown in Fig. [I] It is to be noted that for the pure 
C013 cluster, all the atomic moments are ferromagneti- 
cally coupled and our predicted magnetic moment of the 
ground state structure is in good agreement with the ex- 
perimental value which is 2.0 /is/atom.^ 

In the case of single V atom doped C012V clusters, we 
have again considered the starting geometries of icosa- 
hedral, HCP and cub-octahedral symmetries. However, 
depending upon the position of the singly doped V-atom 
either at the center or on the surface, it can give rise 
to several homotops. Again, the surface atoms of a 13- 
atoms cluster are not equivalent in terms of number of 
neighbors and their relative positions in case of HCP and 
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(0, 0.00 eV) 




(2, 0.17 eV) 



(3, 0.22 eV) 



FIG. 1: (Color online) Structures of optimal HCP, icosahe- 
dron and cub-octahedron of pure C013 cluster (from left to 
right respectively). The first entry in the parenthesis gives 
the isomeric position ( "0" means minimum energy state) and 
the second entry corresponds to relative energy to the mini- 
mum energy state. 



cuboctahedral structures, while for icosahedral structure, 
all the surface atoms are equivalent. This consideration 
will give few more starting geometries. So we have con- 
sidered all these possible geometrical structures and for 
each of them, we have considered the all possible collinear 
spin alignments during relaxation. Upon optimization of 
geometry and spin degrees of freedom, we find that an 
icosahedral structure of total magnetic moment 25 
with V atom doped at the center position, is the MES. 
Now there are several interesting points to note in the 
context of the C012V cluster. First of all, the cohesive 
energy of the MES is considerably higher (by 1.29 eV) 
than that of the pure C013 clusters, indicating its ex- 
ceptional stability. Secondly, the V doped Co-clusters 
are found to prefer the icosahedral growth pattern, in- 
stead of HCP growth pattern, as we observed in case of 
the pure Co-clusters. Therefore, just the single V-atom 
doping in C013 changes the equilibrium structure from 
HCP to icosahedron. This HCP to icosahedron struc- 
tural transition has also been reported previously for Mn 
doping in C013 cluster ™ Finally, the single V-atom likes 
to fit at the center of the icosahedron and it is ferri- 
magnetically coupled with the surface Co-atoms. The 
first isomer is also a central V-atom doped icosahedral 
C012V cluster, with total magnetic moment 23 /i# and 
it is about 0.46 eV above the MES. Center V doped cub- 
octahedral C012V cluster with 23 \±b magnetic moment 
which has energy 0.70 eV higher than the minimum en- 
ergy icosahedral structure, is the second isomer. The 
third isomer is a center V doped icosahedron with total 
magnetic moment of 27 \±b- Center V doping in case of 
HCP cluster is less favorable and appears as the fifth iso- 
mer in our calculation with total magnetic moment 25 /i# 
and energy 0.80 eV above the MES. We have also con- 
sidered the optimal structure of C013 (distorted HCP) 
as the starting structure and substituted the most coor- 
dinated central Co-atom by a V-atom. After relaxation 
considering all possible spin alignments, it is found that 
the shape of the optimal structure remains more or less 
same as that of the optimal HCP C013 cluster and en- 
ergetically, it is the fourth isomer with total magnetic 
moment 23 \±b- The most probable surface V-doping 
has a cub-octahedral structure with total magnetic mo- 



ment 21 but it is about 1.32 eV above the MES. The 
structure of ground state and few isomers for the C012V 
cluster are shown in Fig. [2] 




(0, 0.00 eV) 



(2, 0.70 eV) 



(4, 0.73 eV) 



(5, 0.80 eV) 



FIG. 2: (Color online) Structures of optimal icosahedron, op- 
timal cub-octahedron and two HCP geometries (from left to 
right) for C012V cluster. Green dot represents Co atom, while 
yellow dot represents V-atom. The entries in the parenthesis 
have same meaning as in Fig. [I] 



Because of the icosahedral growth preference of the 
C012V cluster and the presence of the singly doped V- 
atom at the central site, we consider only the icosahedral 
symmetry based geometries as the starting guess for more 
than one V atom doped clusters, in which one V atom is 
always at the center position, while the residual V atoms 
reside on the surface. For C011V2 cluster, the second V 
atom can replace any of the surface Co atoms, as all the 
surface sites of a 13-atom icosahedron, are equivalent. 
After relaxation for all possible spin configurations, we 
find that the MES has total magnetic moment of 19 \±b 
and total cohesive energy of 43.22 eV. The first and sec- 
ond isomers have magnetic moments of 21 /i# and 17 /i# 
and they are 0.10 eV and 0.25 eV above the minimum en- 
ergy state respectively. In the MES, the center V-atom 
is again ferrimagnetically coupled with the surface Co 
atoms and has much lower local magnetic moment as it 
was in the case of the C012V cluster, while the surface V 
atom has the maximum local magnetic moment and it is 
also ferrimagnetically coupled with the other surface Co 
atoms. 




(0, 0.00 eV) 



(5, 0.20 eV) 



(8, 0.31 eV) 



FIG. 3: (Color online) A, B, C represent the three inequiv- 
alent homotops in icosahedral C010V3 structure having one 
V-atom always at the center. The three structures are the 
optimal structures for the respective three types. The entries 
in the parenthesis have same meaning as in Fig. [1] Color 
convention for atoms is same as in Fig. [2] 



For C010V3 cluster, depending upon the different po- 
sitions of the two surface V atoms, there are three possi- 
ble icosahedral structures as shown in Fig. [3j However, 
the structure of type A where the two surface V atoms 
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are closest to each other, appears as the most favorable 
one. This type A icosahedral structure with total mag- 
netic moment of 21 \±b is the minimum energy state with 
ferrimagnetic alignment of the central V-atom and ferro- 
magnetic alignment of the surface V-atoms with the sur- 
face Co-atoms. Both the central and the surface V-atoms 
have small values of local magnetic moments in this case. 
There are several closely spaced isomers of type A icosa- 
hedron with magnetic moments of 19 13 23 hb 
and 17 fi B which are just 0.04 eV, 0.08 eV, 0.09 eV and 
0.15 eV above the MES. We find that the optimal struc- 
tures of type B and type C are 0.20 eV and 0.31 eV above 
the MES, respectively and both of them have same total 
magnetic moment of 13 

For C09V4, there are three V atoms on the surface. 
Considering two surface V atoms at closest to each other 
(following the ground state configuration of C010V3), dif- 
ferent positions of the third surface V atom can give rise 
to the four icosahedral configurations as shown in Fig. 
[4j After optimization, it is found that the type A is the 
most favorable structure where all the three surface V 
atoms are closest to each other and form an octahedron 
with the central V atom. The optimal structure of type 
A icosahedron has total magnetic moment of 15 \±b and 
total cohesive energy of 42.35 eV. Another type A icosa- 
hedrons with magnetic moments of 13 19 21 hb 
and 17 fi B are just 0.03 eV, 0.07 eV, 0.09 eV and 0.13 eV 
above the MES. The optimal type £?, type C and type 
D icosahedrons have total magnetic moments of 15 
13 hb and 15 \±b respectively and they are 0.15 eV, 0.25 
eV and 0.33 eV above the MES. 




(0, 0.00 eV) (6, 0.1 5 eV) (1 1 , 0.25 eV) (1 5, 0.33 eV) 



FIG. 4: (Color online) A, B, C, D represent the four probable 
homotops in icosahedral C09V4 structure having one V-atom 
at the center and two of the remaining three V-atoms are 
closest to each other on the surface, while the position of 
the fourth V-atom is varied. The four structures represented 
above are the optimal structures for the respective four types. 
The entries in the parenthesis have same meaning as in Fig. 
[l] Color convention for atoms is same as in Fig. [2] 

The results of the structural optimizations, therefore, 
can be summarized as, unlike an HCP growth pattern of 
the pure Co cluster, the V doped C013 clusters prefer to 
adopt an icosahedral packing. In such clusters, the most 
coordinated central site is occupied by a V-atom, while 
the residual V-atoms sit on the surface. The surface V- 
atoms like to be closer to each other to form a group, 
thereby imparting more distortion to the structure and 
significantly alter the local surface charge density. The 
central V-site in the MES of all compositions, is always 



ferrimagnetically coupled and has lower magnetic mo- 
ment. On the other hand, the surface V-atoms can be 
ferrimagnetically or ferromagnetically coupled with the 
surface Co-atoms and their local magnetic moments can 
be as low as the central V-site or as high as the sur- 
face Co-site, depending upon the distribution of surface 
charge density in the clusters having different amount of 
V-doping. Our prediction of icosahedral geometry for the 
minimum energy states of V doped Co clusters is in ac- 
cordance with the speculation of a closed shell geometry 
around the central site in Ref. 13 and Ref. 15. 
To have more confidence about the structural trend as 
observed above, we have carried out few ab-initio molec- 
ular dynamics (MD) study for the C012V and C010V3 
clusters. To determine the lowest energy structures, we 
have heated up the clusters to 2000 K (which is close to 
the melting temperature 1768 K for bulk Co and 2183 
K for bulk V). The clusters have then been maintained 
at this temperature for at least 6 ps and then allowed 
to cool again to K. The cooling process has been done 
for 24 ps and maintained a very slow cooling rate of 1 
K per one iteration throughout the process. The main 
results obtained from the MD run are : (a) it prefers to 
adopt the icosahedral pattern with same magnetic align- 
ment of atoms as we predicted from zero temperature 
relaxation, (b) For the C012V cluster, the V atom occu- 
pies the central position, (c) the surface V atoms for the 
C010V3 cluster, prefer to be closest to each other. These 
results provide reassurance and more confidence in our 
zero temperature relaxation. 



IV. STABILITY ANALYSIS 

The observed atomic arrangement of the MES of the 
Coi3_ m V m clusters as described in the previous sec- 
tion depends critically on the balance of several pa- 
rameters like the relative strengths of various kinds of 
bonds present in a cluster structure, the relative atomic 
sizes, the amount of charge transfer between two differ- 
ent species of atoms, the energy gap between the HOMO 
and the LUMO, abundance of states near Fermi energy 
etc. Below, we try to understand the structural stability 
of the clusters in terms of these parameters. 



A. Cohesive energy 

Fig. [5] shows the plot of the cohesive energies of the 
MES of the Coi3_ m V m clusters with increasing V con- 
centration. The plot is with respect to the cohesive en- 
ergy of the MES of C013 cluster. It is seen that the bind- 
ing has increased considerably compared to that of the 
optimal C013 cluster with single V atom doping. How- 
ever, the cohesive energies decrease sharply for the higher 
concentration of V-atoms. For the double and triple V- 
atom doping, the cohesive energies are above the dashed 
line indicating their higher stability compared to the pure 
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FIG. 5: Cohesive energy of the minimum energy structures 
of Coi3- m V m clusters with respect to the cohesive energy of 
minimum energy pure C013 cluster. The dashed line is the ref- 
erence fixed at the cohesive energy of pure C013. Inset shows 
the second difference (A2E) in cohesive energy [as defined in 
Eqn. ([2)] for C012V, Co n V 2 and C010V3. 



C013, while for the fourth V-atom doping, the cohesive 
energy is even lower than that of the pure C013. The 
relative stability among the clusters of nearby concentra- 
tions is more distinct when we plot the second difference 
in cohesive energies in the inset of Fig. [5] We use Eqn. 
[2] to calculate the second difference. A sharp pick in the 
A2E at m = 1 points to the exceptional stability of the 
single V-doped cluster compared to the undoped or more 
than one V-doped clusters. 

In order to understand the stability of the doped clus- 
ters, we have calculated the cohesive energies and the 
bond lengths of the Co-Co, Co-V and V-V dimers as 
shown in Table [TJ It is seen that the V2 dimer is the 
most stable and the bond length of the V2 dimer is also 
about 14% shorter than that of the C02 dimer which 
has the least binding among the three. On the other 
hand, the cohesive energy and bond length of the Co-V 
dimer are intermediate of the C02 and the V2 dimers. For 
the C012V cluster, we have seen that the center doping 
in icosahedral structure is the most favorable. Binary 
clusters are known to have the problem of segregation, 
where the doped atom can segregate at the surface or 
center. To a first approximation, if one of the homonu- 
clear bonds is the strongest, then that species tends to 
be at the center of the cluster P Our finding of central 
V-doping is in accordance with this, as the V2 dimer has 
the strongest binding. The atomic radius of a Co-atom 
and a V-atom, being almost same, the substitution of 
the central Co atom in an icosahedral C013 cluster by a 
V atom leaves the structure almost unperturbed. As for 
evidence, the center to vertex average distance and the 
average distance between two nearby vertices in the opti- 



TABLE I: Cohesive energies and bond lengths of C02, Co-V 
and V2 dimers in the present calculation. For comparison, we 
have also listed the experimental values for C02 dimer (Ref. 
35) and V2 dimer (Ref. 36). 



Dimer 


Cohesive ener^ 


(eV/atom) 


Bond length (A) 




Theory 


Expt. 


Theory Expt. 


Co 2 


1.45 


1.72 


1.96 2.31 


CoV 


1.53 




1.87 


v 2 


1.81 


2.47 


1.72 1.77 



mal icosahedral structure of the C013 cluster are 2.334 A 
and 2.450 A. For the optimal C012V cluster, these values 
are almost same, 2.344 and 2.465 A respectively. There- 
fore the large gain in cohesive energy of the C012V cluster 
over that of the C013 cluster may be due to the enhanced 
cohesive energy of the CoV dimer over that of the C02 
dimer, though an atom in a dimer is quite different from 
an atom at the center of a 13-atom cluster. The V atom 
being at the central position, there are maximum number 
of CoV dimers in V-doped clusters. 

For the clusters doped with more than one V atom, 
the favorable structure is again found to be of an icosa- 
hedral motif, in which the central site is always possessed 
by a V-atom and the other V-atoms lie on the surface. 
The number of the surface V atoms and the V-V bonds, 
therefore, increases with increasing V doping. It favors 
better binding. On the other hand, the presence of the 
V-atoms on the cluster surface distorts the cluster ge- 
ometry and increases the center to vertex average dis- 
tance as well as the vertex-vertex average distance. Ef- 
fectively, the cluster volume increases with the increasing 
V-concentration, which destabilizes the structure. In Ta- 
ble [TTJ we have listed the average distances between cen- 
ter to vertex as well as between two nearby vertices for 
the MES of Coi 3 _ m V m (m = 0-4). We believe that the 
distortion in the structures is due to the charge density 
variation induced by the presence of the V-atoms on the 
cluster surface. The strained cluster structures resulted 
from distortion can be realized from some open bonds in 
Fig. [3] and Fig. [4j Therefore, though the number of the 
V-V bonds increases, which has better binding compared 
to the Co-Co or Co-V bonds, the overall cohesive energy 
decreases monotonically as we go from C012V to C09V4 
by successive doping of V-atoms. It is then obvious that 
the cluster geometry and the distribution of atoms on 
the cluster surface play an important role in deciding the 
cluster stability. 

At this point, it is interesting to study the structure 
and energetics of the pure V13 cluster. Since experiment 
hints symmetric structure for the V13 cluster^, we con- 
sidered an icosahedral geometry and optimized it consid- 
ering all possible spin configurations. Interestingly, the 
decreasing trend of cohesive energies for the Coi3_ m V m 
clusters with increasing V-concentration continues also 
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TABLE II: The average distances in A between center to ver- 
tex and between two nearby vertices for the minimum energy 
structures of all the studied clusters. For C013, the values 
correspond to the optimal icosahedron. 



TABLE III: Positions of HOMO and LUMO in both the spin 
channels and the values of Si and 82 for the optimized struc- 
tures of all compositions. 



Bonds 


C013 


C012V 


C011V2 


C010V3 


C09V4 


center-vertex 


2.334 


2.344 


2.354 


2.375 


2.398 


vertex-vertex 


2.455 


2.465 


2.476 


2.499 


2.509 



at m = 13 i.e pure V13 cluster. The cohesive energy of 
the optimal V13 cluster is found to be lower than that 
of the pure C013 cluster by about 3.73 eV. It indicates 
that V13 would be more reactive than any of the V-doped 
C013 clusters. This tendency may be rationalized by con- 
sidering the case of bulk V and Co, where the bulk V is 
more reactive than Co because of the lower d-band filling, 
relatively higher position of d-band center and larger ri- 
band width of V compared to Co. However, pure V does 
not have importance in catalysis because its high reac- 
tion tendency with absorbate tends to poison the surface, 
without leaving any active site. The desorption energy is 
also high because of strong V-absorbate bond, which is 
again not favorable for catalysis. 



B. Spin gap 

Analogous to HOMO-LUMO gap of a nonmagnetic 
cluster, one can define spin gaps for a magnetic cluster 

as, 



majority 

"homo 

minority 
e HOMO 



minority 
e LUMO 



majority 
e LUMO 



(3) 



and the system is said to be stable if both S\ and 82 are 
positive i.e the LUMO of the majority spin lies above 
the HOMO of the minority spin and vice versa. These 
represent the energy required to move an infinitesimal 
amount of charge from the HOMO of one spin channel 
to the LUMO of the other. So magnitude of spin gaps is 
a measure of chemical activeness of clusters. Higher gap 
means less reactive and vice versa. The positions of the 
HOMO and LUMO in both the spin channels and the 
values of 8\ and 82 for the optimized structures of 13- 
atom V-doped Co clusters of all compositions considered 



here are given in Table III 



It is seen that both the spin gaps Si and 82 are positive 
for all the clusters. Also C012V has the maximum value 
of 5's which again indicates the highest stability of C012V 
cluster compared to the others. Because of this large spin 
gap, C012V has very low reaction tendency towards H 2 
molecules. However, gap values decrease with increasing 
V concentration and therefore reactivity also increases as 
observed experimentally!^ 



Cluster 


Majority channel 


Minority channel 


Spin gap (eV) 




HOMO 


LUMO 


HOMO 


LUMO 


Si 


82 


C013 


-3.60 


-3.35 


-3.48 


-3.34 


0.26 


0.13 


C012V 


-3.93 


-2.61 


-3.39 


-3.34 


0.59 


0.78 


C011V2 


-3.71 


-2.83 


-3.43 


-3.32 


0.39 


0.60 


C010V3 


-3.41 


-2.98 


-3.34 


-3.26 


0.15 


0.36 


C09V4 


-3.46 


-3.02 


-3.33 


-3.17 


0.29 


0.31 



C. Density of states (DOS) 

Fig. [6] shows the spin-polarized total DOS, Co pro- 
jected DOS, V-projected DOS and total d-orbital pro- 
jected DOS for the MES of Coi 2 V, ConV 2 , Coi V 3 and 
C09V4 clusters. It is seen that the total DOS and the to- 
tal d-projected DOS are almost overlapping each other. 
This indicates that the cluster properties are mostly dom- 
inated by the d-electrons which is generally expected for 
transition metal clusters. The trend in structural sta- 
bility and reactivity of the clusters as discussed in the 
previous sections, is also obvious from the nature of the 
total DOS. The majority spin channel of each cluster has 
a gap. This gap is maximum for the C012V cluster and it 
decreases as the number of surface V-atoms increases. In 
the minority spin channel, there is finite amount of states 
at the Fermi energy and these are contributed solely by 
the exterior surface atoms, as the central atom does not 
have any contribution at the Fermi energy (cf. Fig. [7]). 
It is also seen that the states are more localized in case 
of the C012V cluster, while they are gradually spread- 
ing out and therefore the DOS height, specially in the 
majority spin channel, decreases with the increasing V- 
concentration. At this point it is relevant to mention 
a common feature observed in case of extended surface. 
The gaseous molecules absorb well on clean surface of 
an early transition metal as well as on the surface with a 
multilayer of late transition metal, but not on the surface 
with a monolayer of late transition metal. For example, 
the photoemission experiments by El-Batanouny et aP^ 
showed that a H atom adsorbs well both on the clean 
Nb(110) surface as well as on the surface with a mul- 
tilayer of Pd, but does not adsorb on Nb(110) surface 
with a monolayer of Pd and this is due to the decrease 
of density of states of ^-electrons of Pd near the Fermi 
level. Similar behavior had been observed also in case of a 
monolayer of Pd on W(110) surface, in the experiment of 
CO-adsorptionP^ Although, a cluster of 13 atoms is sub- 
stantially different in its property from that of the surface 
or bulk, one may resemble the stability of the C012 V clus- 
ter towards H-adsorption with that of extended surface, 
where the central V atom of C012V corresponds to an 
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FIG. 6: (Color online) Total DOS per atom, Co projected DOS per atom, V projected DOS per atom and total d-projected 
DOS per atom for the optimal structures of C012V, C011V2, C010V3 and C09V4. The smearing width is fixed at 0.1 eV. Vertical 
line through zero is the Fermi energy (Ef). 



early transition metal substrate and the exterior twelve 
Co atoms correspond to a late transition-metal mono- 
layer. 



In order to see the chemical activity of the central V- 
atom, we have plotted in Fig. [7j the projected total DOS 
of the central V-atom for each of the Coi3_ m V m clusters. 
First of all, it is seen that there is no finite states at 
the Fermi energy, meaning that the central V-atom is 
not chemically active. Also each of the majority and 
minority spin channels has two large pecks: one is deep 
below the Fermi energy and another is above the Fermi 
energy. However, the peck heights gradually decrease, 
the states are broadened out and shifted towards higher 
energy with increasing V-concentrations. It is therefore 
indicating that the presence of surface V-atoms induces 
some sort of chemical activeness to the central atom. 



V. CHEMISORPTION WITH H 2 MOLECULES 



In order to gain some understanding about the cluster 
chemical reactivity, we have investigated the chemisorbed 
structures of the Coi3_ m V m clusters upon H2 adsorb- 
tion. To check the robustness of our chemisorption cal- 
culations involving H-atoms, we have first calculated the 
cohesive energy, bond length and vibrational frequency 
of H2 dimer. Our calculated values have been listed in 
Table [IV| These values are typical for gradient-corrected 
calculations of H2, which have been done before^ and 
they agreed reasonably well with experiment.^ 

We have then performed an exhaustive search for the 
MES, taking H2 at different possible places on the MES of 
the corresponding bare cluster for each composition. Fig. 
|8]shows our calculated lowest energy chemisorption struc- 
tures after H 2 absorbsion for clusters of all compositions. 
First thing to notice is that H2 molecule chemisorbs dis- 
sociatively in each case, i.e the distance between the two 
H-atoms in the chemisorption structures is much larger 
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TABLE IV: Theoretical values of cohesive energy, bond length 
and vibrational frequency for H2 dimer in the present calcu- 
lation. Experimental values in Ref. 41 are also given for 
comparison. 



Cohesive Bond length (A) Vibrational 

energy (eV) frequency (cm -1 ) 
Theory 4.520 0.752 4339 
Experiment 4.750 0.741 4395 



than the H-H bond length of an isolated H2 molecule. 
The chemisorption gives rise to moderate perturbation 
to the structures and makes them more symmetric (i.e 
surface bonds are now less dispersive) compared to the 
parent structure without hydrogen. There are three pos- 
sible ways for the H-atoms to be adsorbed on each cluster 
: on top of an atom (one fold), at a bridge position be- 
tween two atoms (two fold) and at the hollow site of a 
triangular plane on the cluster (three fold). Again, as the 
surface contains two species of atoms for C011V2, C010V3 
and C09V4, then the one fold position can be on top of a 
surface V-atom or on top of a surface Co-atom. Two fold 
position can be the top of a V-V bond, Co-V bond or a 
Co- Co bond which can be nearer to or away from the V- 
site. Similarly, in triangular plane, there are several pos- 
sibilities : (i) all the three atoms can be V-atoms (only 
possible for C09V4), (ii) one Co atom and two V-atoms, 



(Hi) two Co-atoms and one V-atom, and (iv) all three are 
Co- atoms. We have considered all these possible combi- 
nations during optimization. It is, however, seen that 
in each of the optimized structures, H-atoms absorb at 
the hollow site on the surface and for m > 2, they pre- 
fer the association of the local V atom. For example in 
the C09V4H2 cluster, one H-atom absorbs at the hollow 
site of V-V-V triangular plane and the other H atom on 
top of a V-V-Co triangular plane. For C011V2H2, the 
H-atoms appear to absorb at the bridge positions, but 
they are inclined with an angle such that the absorption 
tends toward a three fold configuration. The preference 
of more coordinated hollow site is likely due to geomet- 
ric arrangements. It allows the H-atom to interact more 
with the V or Co atoms. On the other hand, the V-site 
preference of H-atom is due to the formation of stronger 
s-d bond with V-atom compared to that with Co-atom. 

The higher reactivity of the C013 cluster compared to 
the C012V cluster, while both of them have the same 
twelve Co- atoms on the surface, can be understood from 
charge transfer analysis. We have computed the charge 
enclosed within a sphere around each surface atom of 
the optimal structures of C013H2 and C012VH2 clus- 
ters. It is seen that the amount of charge on the sur- 
face Co-atoms in the C013H2 cluster is larger by about 
0.2 (in unit of electronic charge) than that of C012VH2 
and it is mainly contributed by the d-electrons of surface 
Co atoms. Therefore stronger 3d-ls interaction between 
surface Co-atoms and chemisorbed H-atoms in case of 
C013H2 increases its reactivity. This result is in accor- 
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FIG. 8: (Color online) The calculated minimum energy chemisorption structures with H2 on the minimum energy structures of 
Coi3- m V m (m = 0-4). It is clearly seen that hollow site on the surface is preferred by chemisorbed hydrogen. Color convention 
is same as in Fig. [2] with pink colored dots representing H atom. 
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FIG. 9: The calculated chemisorption energy (using Eqn. [4] 
for Coi3-mV m H 2 clusters. 



dance with that of Fujima et al 42 who also predicted the 
stronger interaction between Is of H atom and the 3d 
orbital of the surface Co atom in the C013H2 cluster 
compared to that in the C012VH2 cluster using density 
of state analysis. They showed that the anti-bonding 
orbital component between the H-ls electron and the 
3d electron of the surface Co-atoms shifts further away 
above the HOMO in case of C013H2 compared to that 
of C012VH2, while the bonding orbital components for 
both the clusters stay almost at the same energy posi- 
tion below the HOMO. Extending the charge transfer 
analysis for the Coi3_ m V m H2 clusters with m = 2,3 and 
4, we find that the average charge per surface V-atom is 
gradually increasing with the increase of V-concentration 
which indicates the formation of stronger V-H bonds at 
the surface with increasing V-doping. Average charge 
per surface Co-atom for these three clusters, however, 
remains almost the same. 

Fig. [9] shows the plot of chemisorption energy with 
increasing V-concentration, where the chemisorption en- 



ergy is defined as 

D e (E) = £(Co 13 _ m V m ) + E(R 2 ) - £(Co 13 - m V m H 2 ) 

(4) 

It is a positive quantity and gives a measure of binding 
strength of H2 molecule with the cluster. The plot shows 
a minimum for C012VH2 again indicating the lowest 
binding efficiency of C012V with H. However, chemisorp- 
tion energy increases with increasing V concentration. 
The source of this chemisorption energy is the cluster re- 
arrangement energy (ie the energy change due to the ge- 
ometrical rearrangement of the cluster upon chemisorp- 
tion) and the efficient cluster-absorbate bonding in pres- 
ence of V due to more efficient s-d hybridization. 

VI. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the geometric and elec- 
tronic structures of V doped C013 clusters and their 
chemisorption towards hydrogen molecules using first 
principles density functional calculation. The lowest en- 
ergy structures of all compositions prefer to have icosa- 
hedral geometry, unlike hexagonal symmetry preference 
of the pure Co clusters. For the C012V cluster, the single 
V atom prefers to sit at the central site, thereby guarded 
by all the surface Co atoms and cannot participate di- 
rectly in the chemisorption reaction. Consequently reac- 
tivity of C012V becomes very less. On the other hand, for 
more than one V atom doping, the additional V-atoms 
reside on the surface and come in direct contact with 
chemisorbed H-atom and reactivity increases. Our calcu- 
lated spin gaps, density of states and charge transfer anal- 
ysis explain nicely the stability of clusters and their ten- 
dency towards chemisorption. In the chemisorbed struc- 
tures, H- atoms adsorb dissociatively at the more coordi- 
nated hollow sites and they prefer V-site association due 
to stronger 3d-ls hybridization. To have better insight 
into the chemisorption reaction, one needs to study the 
transition states for the optimal cluster of each composi- 
tion and we believe that calculation of activation barriers 
will also lead to same conclusion as we have predicted 
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here. 
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